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SAMPLING METHOD FOR A SPREAD SPECTRUM COMMUNICATION SYSTEM 

The present invention relates to a method for decoding signals 
sent over a bandwidth-expanding communication system, said method 
5 allowing to perform both channel estimation and signal detection. The 
present invention also relates to a receiver using said method as well as to a 
communication system including said receiver. 

A communication system is defined to be bandwidth-expanding 
if the transmitted signal occupies a bandwidth much larger than the 

10 minimum bandwidth required for sending the information contained in 
said transmitted signal. Bandwidth-expanding communication systems, such 
as Code Division Multiple Access (CDMA) communication systems or Ultra- 
Wide Band (UWB) communication systems for example, possess many 
intrinsic advantages over conventional communication systems, among 

15 which are selective addressing capability, high-resolution ranging, low- 
density power spectra for signal hiding, interference rejection, etc. 

In a CDMA communication system, for example, bandwidth 
expansion is accomplished by means of a spreading code, often called a 
coding sequence, which is independent of the information data to be sent. 

20 Each symbol (information bit) of a signal sent over a CDMA communication 
system is coded with a coding sequence. The number of chips Nc in the 
coding sequence is referred to as the spreading factor or the processing 
gain of the communication system. The chip rate of the transmitted signal 
is thus Nc times higher than its symbol rate. Coding an information signal 

25 by the coding sequence therefore expands the signal's bandwidth by the 
value of the spreading factor Nc. 

The spreading factor of a bandwidth-expanding communication 
system is typically a large number so as to provide said system with the 
above mentioned advantages over conventional communication systems. 
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In a CDMA communication system for example, the spreading 
factor Nc, or processing gain, determines the maximum number of users of 
the system. Namely, each user of the communication system is assigned a 
unique coding sequence which must have a low cross correlation with the 
5 other users' coding sequences, in order to minimize interference between 
users. In the ideal case, the user's codes are chosen from a set of orthogonal 
sequences. Therefore, the more users in the communication system, the 
longer these sequences must be, thus increasing the spreading factor No 

As shown in Figure 1, the signals Xk(t) sent by the users k over the 

10 bandwidth-expanding communication system are first encoded with a user 
specific coding sequence Sk(t) by encoders 50 and then transmitted over a 
transmission channel c comprising emitters 60 and at least one receiver 61 . 
The emitters 60 and the at least one receiver 61 can be for instance radio 
emitters and receivers working with amplitude or frequency modulation of 

15 a carrier signal. The received wideband signal y(t), which Is a superposition 
of all signals 10 transmitted by the active users k over the transmission 
channel c must be decoded at the receiver's side in order to retrieve the 
signal Xk(t) sent by each user k of the system. For example, a conventional 
analog CDMA receiver 70 comprises a bank of analog matched filters fk, 

20 each matched to a specific user's coding sequence Sk(t). The received 

wideband signal y(t) is filtered in parallel by each filter fk. The information 
bits of the signal Xk(t) sent by each user k can then be retrieved from the 
output signal of the corresponding filter fk. This works well if there are few 
users with almost orthogonal coding sequences and if all transmitted 

25 signals have equal power levels. However, in the case when these power 
levels are not equal, the standard detector becomes almost useless. This Is 
known as the near-far problem, and is currently usually alleviated by using 
for example complex power-control schemes, multiuser detection schemes, 
RAKE receivers, adaptive antenna arrays or 2-D RAKE receivers. 

30 Current trend in design of communication systems is to 

implement as many parts of receivers as possible in digital hardware. The 
main reason is the fact that digital circuits are cheaper, faster and smaller. 
These receivers however imply that the received continuous-time signal 
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must be first converted into a discrete-time signal by sampling. Defining 
the bandwidth of the received signal to be between the first two zeros in 
the signal's frequency spectrum, then according to the Shannon's sampling 
theorem, the sampling frequency must be at least twice the maximum 
5 frequency of the signal. The required minimal sampling frequency defined 
by this theorem, commonly referred to as the Nyquist frequency, increases 
with the bandwidth of the signal to be sampled. 

For example, in a CDMA communication system, this means that 
the sampling frequency must be greater or equal to twice the chip rate of 

10 the received signal. In practice, however, in order to perform good signal 
detection even in noisy or fading multipath environments, the sampling 
frequency used by most current receivers is at least four times higher than 
the chip rate, thus requiring said receivers to include very fast and 
therefore complex and expensive analog sampling devices. Faster sampling 

15 devices have furthermore high power consumption and can hardly be 
integrated in mobile battery powered receivers. 

In the case of a CDMA communication system, as already 
mentioned, the near-far problem can be alleviated by using for example 
multiuser detection schemes, RAKE receivers, adaptive antenna arrays or 2- 

20 D RAKE receivers, thus improving the system's performances. These 

methods can be implemented in digital receivers as well. However, they 
include sophisticated signal processing techniques that require exact 
knowledge of one or several parameters of the transmission channel, such 
as relative delays of different users' signals with respect to a reference in 

25 the receiver, amplitude attenuations of these signals along the different 
propagation paths, their direction-of-arrival at the receiver, etc. 

Estimation of these parameters, referred to as channel 
estimation, is crucial in all bandwidth-expanding, or spread-spectrum 
communication systems, but at the same time is the hardest part of system 
30 design. In radio communication systems, for example, channel estimation, 
which can be performed during a dedicated training phase, represents the 
major problem, particularly in the case when a mobile system is subject to 
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multipath fading, that is, when the transmitted signal follows different 
propagation paths from the emitters to the receiver, each path having 
different attenuation characteristics, delays and directions-of-arrival. The 
existing solutions typically require sampling the received signal at a high 
5 sampling frequency at the receiver's side, the sampling frequency being 
usually a multiple of the signal's chip rate, thus generating an excessive set 
of samples during the training phase. The propagation parameters are then 
estimated from this set of samples, by solving an optimization problem that 
involves large dimension matrices and complicated temporal and/or spatial 
10 operations, thus making this approach computationally intensive and 
sometimes unaffordable in real-time equipment. Once the propagation 
parameters have been estimated for the signals transmitted by all users 
over the transmission channel, the digital receiver can proceed with signal 
detection, that is, retrieving the information sent by each user. 

15 An aim of the present invention is to provide a method for 

decoding a signal sent over a bandwidth-expanding communication system, 
that is, a method for performing both channel estimation and signal 
detection, said method allowing the design of digital receivers comprising a 
sampling device with reduced complexity and cost compared to the 

20 complexity and cost of the sampling devices of current digital receivers, 
while retaining equivalent decoding performances. 

Another aim of the present invention is to provide a decoding 
method for signals sent over a bandwidth-expanding communication 
system, allowing the design of digital receivers with reduced computational 
25 requirements compared to the computational requirements of current 
digital receivers, while retaining equivalent decoding performances. 

These aims are achieved by a decoding method and a related 
receiver comprising the features of the corresponding independent claims, 
features of preferred embodiments being described in the dependant 
30 claims. 
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In particular, the method according to the invention allows the 
decoding, that is, both channel estimation and signal detection, of a signal 
sent over a bandwidth-expanding communication system from a set of 
sampled values generated by sampling the received wideband signal with a 
5 sampling frequency lower than the sampling frequency given by the 

Shannon's sampling theorem, but greater than the rate of innovation p of 
said sent signal, thus allowing a receiver using this method to sample the 
received signal at a sub-Nyquist rate and still allowing an exact 
reconstruction of the sent signal. 

10 With the method according to the invention, all the necessary 

computational steps for channel estimation and/or for signal detection are 
performed on the basis of a significantly reduced set of sampled values 
compared to the set of sampled values needed by current digital receivers, 
thus allowing a significant reduction of the computational requirements of 

15 the digital receiver using the method according to the invention compared 
to the computational requirements of current receivers having equivalent 
performances. 

In particular, during the training phase, the method according to 
the invention allows estimation of the propagation parameters, that is, 

20 relative delays of different users' signals with respect to a reference in the 
receiver, propagation coefficients along the different propagation paths in 
the case when the communication system is subject to multipath fading 
and direction-of-arrival of the different signals, from a significantly reduced 
set of sampled values compared to the set of sampled values needed by 

25 current digital receivers. 

The minimal sampling frequency required by the method 
according to the invention is determined by the rate of innovation p of the 
received signal. The rate of innovation p of a signal is defined as the 
number of degrees of freedom of the signal per unit of time. For example, 
30 in the case of a periodic stream of Z weighted pulses (Diracs) per period T, 
the signal is fully specified by the Z amplitudes and the Z locations of the 
Diracs of one period T. That is, the number of degree of freedom of this 
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signal is 2Z per period T, and its rate of innovation is therefore p=2Z/T. In 
the case of a signal sent by a user k over a bandwidth-expanding 
communication channel, given a known spreading function, the wideband 
signal transmitted by said user k is entirely determined by the sent 
5 information or sent symbols. In the single-user case, the number of degrees 
of freedom of the received wideband signal is therefore equal to the 
number of degrees of freedom of the sent signal, that is the number of 
symbols sent, and the rate of innovation p of the received signal is equal to 
the number of symbols sent by the user k per unit of time. The rate of 

10 innovation p of the received signal y(t) is therefore equal to its symbol rate 
1/Tb, which, as explained before, is lower than its chip rate 1/Tc by the 
spreading factor. In the multi-user case, the rate of innovation of the 
received signal is equivalent to the addition of the rates of innovation of 
the K signals sent by the K users k. The rate of innovation p of the received 

1 5 signal y(t) is then equal to the number K of symbols sent by the K users k 
per unit of time, the rate of innovation p of the received signal y(t) is 
therefore equal to its information rate K/Tb. The number K of users k being 
typically significantly smaller than the spreading factor Nc, the rate of 
innovation p of the received signal is still significantly lower than its chip 

20 rate 1/Tc. 

The invention will be better understood with the help of the 
figures 1 to 5. 

Figure 1 which has been previously described illustrates a CDMA 
system including a known analog receiver. 

25 Figure 2 diagrammatically represents a bandwidth-expanding 

communication channel, for example a CDMA system. 

Figure 3 illustrates the concept of multipath propagation from 
the emitter to the receiver within a transmission channel. 

Figure 4 represents the steps performed by the processing device 
30 w of a digital receiver r using the method according to the invention. 
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Figure 5 illustrates the principle of an array of antennas. 

In the following of the specification and claims, the terms 
communication system and communication channel are equally used, unless 
otherwise specified. 

5 In its preferred embodiment, the decoding method according to 

the invention is applied to the decoding of a signal sent over a CDMA 
communication channel, said channel including a receiver using the 
method according to the invention. 

Figure 2 diagrammatical ly represents a CDMA system with a 

10 number K of users k operating over a transmission channel c. The signal x(t) 
sent by each user k is encoded with the user specific coding sequence Sk(t) 
by an encoder 50, thus generating a wideband signal 10. This wideband 
signal 10 is transmitted to the receiver r, for example using frequency or 
amplitude modulation of a carrier signal. At the receiver r, the received 

15 signal y(t) after demodulation is composed of the superposition of all 

wideband signals 10 transmitted by the K users k. The signal is filtered with 
a lowpass filter f, that is, it is convolved with a filtering signal (p(t), for 
instance a bandlimited sine signal. The filtered signal 12 is then sampled at 
a regular sampling frequency fs=^fTs, resulting in the time-continuous 

20 sampled signal ys(t). The sampled values y[nTs] are extracted from this time- 
continuous sampled signal ys(t) by a continuous-to-discrete converter d and 
the spectral values Y[m] of the time-continuous signal ys(t) are computed 
from these sampled values y[nTs], for instance by using the Fast Fourier 
Transform (FFT) method. The spectral values Y[m] are then used by the 

25 processing device w either to estimate the values of the propagation 
parameters of the transmission channel c during the training phase, or to 
retrieve the information sent by each user k in the detection phase, as 
explained below. 

Figure 3 diagrammatically represents the functioning of the 
30 processing device w, first during the training phase 30 and then during the 
subsequent detection phase 40. 
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During the training phase 30, the processing device w solves a 
one-dimensional estimation problem including a set of spectral values Y[m] 
of the sampled signal ys(t), the known spectral values Sk[ml of each user's 
coding sequence Sk(t) and the known values of the symbols bk^^^of the 
5 training sequences bkt sent by each user k, in order to determine the 
propagation parameters of the transmission channel c. In the preferred 
embodiment of the method, the estimated propagation parameters are the 
delays xk^'^ and the amplitude attenuations ak^'^ induced by the transmission 
channel c on the signals sent by the users k, which in this phase are the 

10 coded training sequences bkt, and transmitted along the different 
propagation paths I. The spectral values Sk[m] of each user's coding 
sequence Sk{t) can be computed by sampling said coding sequence Sk(t) at 
the same sampling frequency fs as the one used for sampling the received 
signal y(t), and applying, for instance, the FFT method to the generated set 

15 of sampled values Sk[nTs]. These spectral values Sk[nTs] can be stored for 
instance in a memory in the receiver r 

During the detection phase 40, the processing device w solves a 
second estimation problem including a second set of spectral values Y[m] of 
the sampled signal ys(t), the previously calculated values of the propagation 
20 parameters ak^'^ and xk^'^ and the known spectral values Sk[m] of each user's 
coding sequence Sk(t), in order to determine the information sent by the 
users k, that is the value of the sent symbols bK. 

The concept of multipath fading transmission channel is 
illustrated by figure 4 which shows how a signal transmitted by a user k can 

25 follow different propagation paths I with different attenuation properties 
and propagation delays. Each copy of the transmitted signal 10 reaching 
the receiver r through a different propagation path I has therefore a 
possibly different complex amplitude attenuation ak^*^ and a possibly 
different relative delay xk^'^ with respect to a reference in the receiver r. The 

30 estimated maximal number of paths I that each transmitted signal 10 can 
follow within a transmission channel is at most equal to L The transmission 
channel is assumed to vary slowly, which means that its propagation 
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parameters ak^'^ and xk^'^ are considered constant during the whole training 
phase. 

In the following, the preferred embodiment of the method 
according to the invention will be described in more detail and it will be 
5 demonstrated how the method according to the invention allows the 
decoding of a signal sent over a bandwidth-expanding communication 
system, by taking samples of the received signal y(t) with a sampling 
frequency fs lower than the sampling frequency given by Shannon's 
sampling theorem, but greater than the rate of innovation p of said signal. 

10 With the signature sequences Sk(t) (si(t), S2(t), SkW) assigned to 

each of the K users k and with bi, bz, bK the symbols bk sent by the K 
users k during the duration Tb of a symbol bk, and assuming that there are 
at most L propagation paths I for each signal transmitted by a user k, the 
received signal y(t) during this period of time, after demodulation, can be 

15 expressed as: 

y(t)=bila/^>si(t-t/^V ai^^^Si(t-Ti^^V...+ a/^^Si(t-T/^^)] 
+ b2[a2'%it^T2^'')^ a2^2>S2(t-T2^^>)+...+ a^^^^s^ft-x^^^^)] 
+ ... 

+ bK[aK<^^SK(t-XK^^V aK^'^SK(t-XK^2V...+ aK^^>SK(t-XK^'^)] 

20 that is, the received signal y(t) is a superposition of multiple 

copies of the signals transmitted by the K users k over the transmission 
channel c, each copy having a possibly different amplitude attenuation ak^'^ 
and a possibly different relative delay xk^'^ with respect to a reference at the 
receiver r. The received signal y(t) can be written more compactly as 

25 y(t) = tb,ta^.(t-x^") 

k=1 1=1 

where tk^'^ denotes the delay with respect to a reference at the 
receiver r of the signal transmitted by user k along the propagation path 1, 
and ak^" denotes a complex amplitude attenuation induced by the 
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propagation path I that includes contributions from the attenuation of the 
transmission channel c and phase offset. 

At the receiver r, the received signal y(t) is filtered with a lowpass 
filter f. In the preferred embodiment of the method of the invention, the 
5 received signal y(t) is filtered with the filter f having the impulse response 
cp{t), where (p(t) is for instance a sine signal of bandwidth j-Mcoor McooL where 
coo=27r/Tb. 

The sampling period T, is chosen such that Tb/Ts > 2M+1. A set of 
at least 2M+1 samples y[nTsl is thus generated during each symbol period 
10 Tb. From the set of samples, 2M+1 spectral values Y[m], where m€[-M,Ml, 
can be computed, for instance with the use of the Fast Fourier Transform 
(FFT) method. 



The delays xk^'^ of the transmitted signals along the different 
15 propagation paths I appear in the spectral coefficients Y[m] as phase 
delays, while the complex amplitude attenuations ak^'^ appear as the 
weighting coefficients. 




k=1 



1=1 



The above equation can be written more compactly as 



K 



mk 



k=l 



where Cmk are given by 




or in a matrix form 
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= C b = Y 



During the training phase 30, the values of the symbols bk sent by 
the users k are known and the linear system above has to be solved for the 
elements of the matrix C, in order to estimate the delays xk^'^ and amplitude 
attenuations ak^'^ of the transmitted signals 10. 



However, the linear system above alone, generated from the 
sampled values y[nTs] of the received signal y(t) taken during one symbol 
duration Tb, is not sufficient to solve for all the elements of the matrix C 
Since the matrix C is of a size M x K, it is necessary to have at least K such 
10 equations, that can be obtained by considering at least K set of samples 
y[nTs] taken from the received signal y(t) during at least K symbol durations 
Tb. In other words, in order to solve for C given b and Y, each user k has to 
send a training sequence bkt, where 



K 

bat 




b<'> 


b<^» 
bf 


- b<;' 

... b<"'; 













15 The length H of each training sequence bkt, which is the number 

H of symbols bk^^^ is at least K, where K denotes the number of users k. 

The equation system to be solved by the processing device w of 
the receiver r can thus be written as 



20 



-21 



-12 



-22 



^2K 



b<" 



(2) 



bl 
b<^' 



AP? b<« 



yd) [2] 



Y»'[11 
Y(2){2] 



Y*«[2] 



sC B = 



wo 02/(178204 



PCT/EP(>2/03390 



This can be solved for C as 
C = Y 

The above equation provides a unique solution as long as B is a full rank 
matrix. B is the matrix made up of the known training sequences bkt of all 
5 users k and, as long as these training sequences bkt are linearly 

independent, the matrix B will be a full rank matrix. In other words, each 
user k must send a determined training sequence bkt, stored for instance in 
a memory storage device in the corresponding encoder 50, linearly 
independent from any training sequence bkt sent by any other user k. 

10 From the above equation all the coefficients Cmk can be 

calculated. 



Since the spectral values 5k[m] of the coding sequences Sk(t) are 
known weighting coefficients, a new matrix D can be defined by 




Cn/S,[1] 

c,,/S,[2] 



C22/SJ2] 



Cik/Sk[11 
C,,/SJ2] 



CM2/SJII 



1; 



with dmk =Cmk/Sk[ml. The coefficients dmk are thus given by 




1=1 



With the number M>2L then the k-th column of D provides 
sufficient information to solve uniquely for the channel parameters of the 
20 user k. For example, considering the first column of D: 
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Only the unknown propagation parameters of the first user 
appear in the above set of equations. In other words, the propagation 
parameters ik^'^ and ak^'^ of user k can be perfectly estimated, in the 
5 noiseless case, from the k-th column of matrix D. The problem of multiuser 
channel estimation is therefore decomposed into a series of one- 
dimensional estimation problems that can be efficiently solved for instance 
using known 1-D subspace methods for harmonic retrieval. 

It is thus demonstrated that the method of the invention allows, 
10 during the training phase 30, the estimation of the propagation 

parameters xk^'^ and ak^'^ for all users k from a lowpass filtered version of the 
received wideband signal y(t) sampled at a sub-Nyquist rate, which implies 
that the sampling frequency fs as well as the computational requirements 
at the receiver's side can be significantly reduced in comparison to the 
15 sampling frequency and the computational requirements induced by 
current methods. 

The number of samples M needed per symbol bk*^^ depends only 
on the number L of multipaths I, while the number H of symbols bk^^^ in 
each training sequence bkt depends on the number K of users k that need 
20 to be synchronized. In the noiseless case, the method according to the 
preferred embodiment of the invention leads to a perfect estimation of 
channel parameters ak^'^ and tk^'^ of all users k by taking only M=2L+1 
samples Y[m] per symbol bk^^^. This means that the sampling frequency fs 
must be at least greater than 2L times the symbol rate: 



25 



wo 02/U78204 



PCT/EP«)2AU390 



14 

During the training phase 30, the values of the sent bits bk^' of 
the training sequences bkt being known, the degrees of freedom of the 
received signal y(t) are the set of unknown propagation parameters, that is 
the KL delays ik^'^ and the KL amplitude attenuations ak^'^ that the K 

5 transmitted signals 10 can take along their L propagation paths I. The 
number of degrees of freedom of the received signal y(t) is thus 2KL. Since 
the propagation parameters are constant over the entire training phase 30, 
and since at least K bits bk^^^ of the training sequences bkt are known, the 
number of degrees of freedom of the received signal y(t) stays constant 

10 over the entire training phase, which is at least KTb long. The rate of 
innovation p of the received signal y{t) can therefore be expressed as 

^ 2KL ^ 2L 
^'KT, ^T, 

Thus it is demonstrated that the sampling frequency has to be 
higher than the rate of innovation of the signal. 

15 _ = f,>p = _ 

In a CDMA communication system, the symbol rate 1/Tb is 
typically considerably lower than the received signal's chip rate l/To The 
inventive method thus significantly reduces the required sampling 
frequency fs in comparison to current methods typically requiring digital 
20 receivers to sample the received signal at least at its chip rate 1/Tc. 

Consequently, the number of sampled values required by the 
inventive method for estimating the propagation parameters is also 
significantly reduced, reducing the complexity of the estimation problem to 
be solved and thus reducing the computational requirements for the 
25 processing device w of the digital receiver r compared to the computational 
requirements for receivers using prior art methods. 



wo 02/078204 



PCT/EP02/03390 



15 

As already mentioned, the training sequences bkt sent by the 
users k during the training phase 30 have to be carefully chosen so as to 
ensure the good performance of the inventive method for channel 
estimation. Therefore, the training sequence bkt to be sent by each user k is 

5 chosen from a set of training sequences specifically designed to optimize 
the performance of the channel estimation and stored in the corresponding 
encoder 50. Preferably, the encoder 50 can store at least two different 
training sequences bkt from this set During the training phase 30, the user 
k sends the stored training sequence bkt most adapted to the current 

10 characteristics of the communication channel, basing for instance on a 
feedback message from the receiver r. 

In the noisy case, in order to obtain a more precise estimation of 
the propagation parameters ak^'^ and tk^'^ the sampling frequency fs needs 
to be increased, while the length H of the training sequences bkt remains 
15 the same. In the presence of noise, the sampling frequency fs required for 
good estimation accuracy depends on signal-to-noise ratio (SNR), yet in 
most cases encountered in practice that sampling frequency fs is still far 
below the chip rate 1/Tc. 

In a variant preferred embodiment, the method described above 
20 during the training phase 30 is applied to an asynchronous CDMA 
communication system where the receiver r has an arbitrary timing 
reference that is not aligned to transmitted symbol bk^^^ boundaries. 
Therefore, in the general case, the sample values ys[nTs] from which the 
Fourier series coefficients Y[m] are computed, will actually contain the 
25 information about two consecutive bits in the training sequence, provided 
that a maximum delay spread of each user k is less than a symbol duration 
Tb- As a result, the computed Fourier series coefficients Y[m] may not 
correspond to the actual values, in particular in the case when the adjacent 
bits have different signs. It is thus desirable to repeat every bit bk^^^ in the 
30 training sequence bkt twice during the training phase, to make sure that 
the Fourier series coefficients Y[ml are computed accurately. This does not 
imply that the derived equations should be modified. All the equations 
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remain the same, the only difference is that every bit bk^*^^ in the training 
sequences bkt is sent twice. 

Once the propagation parameters ak^'^ and ik^'^ of the transmission 
channel c for all users k have been estimated, the processing device w of 
5 the receiver r proceeds to the detection phase 40, that is, the detection of 
the signal x(t), or the symbols bk, sent by each user k. 

According to the preferred embodiment of the invention, the 
detection phase 40 is performed using known multiuser detection methods 
such as Minimum Mean Square Error (MMSE) or decorrelating detector. The 

10 received signal y(t) however is, as in the training phase 30, first filtered with 
a lowpass filter f and then sampled at a sub-Nyquist rate, that is, at a 
sampling frequency fs lower than the chip rate 1/Tc of the received signal 
y(t), but higher than its rate of innovation p. The inventive method base on 
the surprising finding that a set of sampled values generated from a signal 

15 at a sub-Nyquist rate is sufficient to reconstruct or to decode said signal, if 
said set of sampled values is generated at a sampling frequency fs higher 
than the signal's rate of innovation p. 

During the detection phase 40, the receiver r receives a second 
signal y(t) made of the superposition of multiple copies of information 

20 signals x(t) sent over the communication channel by the users k, each copy 
of each information signal x(t) having a previously estimated amplitude 
attenuation ak^'^ and a previously estimated relative delay xk^'^ with respect 
to a reference at the receiver r. The information signals x(t) sent by the 
users k are made of symbols bk. These symbols bk are then coded with the 

25 respective coding sequence Sk{t), and transmitted to the receiver r over the 
transmission channel c. The second signal y(t) is first filtered and then 
sampled at a sampling frequency fs lower than the sampling frequency 
given by the Shannon theorem, but higher than rate of innovation p of the 
second signal y(t), which is, as described above in the multi-user case, equal 

30 to the signal's information rate K/Tb, thus generating a second set of 
sampled values y[nTs]. Signal detection is then performed with known 
detection methods such as for instance Minimum Mean Square Error 
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(MMSE) or decorrelating detector, using the values of the second set 
sampled values y[nTs] and the previously estimated propagation parameters 
ak<'^ and t^^'^ 

In a variant preferred embodiment of the method, the training 
5 phase 30 is initiated as described previously. The coefficients Cmk of the 
matrix C are computed with 

C = Y 

However, once the matrix C is computed, the system proceeds to 
the subsequent detection phase 40 without solving the series of one- 
10 dimensional estimation problems for estimating each propagation 

parameter ak^'^ or z^^^K During the detection phase 40, the second signal y(t) 
is sampled as described in the previous embodiment of the method, thus 
generating a second set of sampled values y[nTs] from which a second set of 
spectral values Y[ml can be computed with for instance the FFT method. 



15 



The linear system 



C2I <^22 
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b = Y 



is then solved for b in order to determine the value of the 
symbols bk, with all coefficients Cmk being known from the synchronization 
phase. Therefore, the vector containing the bits sent by each user can be 
20 computed as 
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In other words, it is possible to detect user's signals without 
knowing specific users' signature sequences and without even estimating 
the propagation coefficients aic^'* and xk^^\ The receiver only has to estimate 
the matrix C during the training phase, while during the detection phase 
5 the vector of bits b is simply given by the linear system above. This 
preferred embodiment of the method thus significantly reduces the 
computational requirements for the computing device w of a receiver r 
using said method. However, it can only be satisfactorily applied to quasi- 
synchronous CDMA communication systems where the timings of the users 
10 k are a priori known to the receiver, and where the maximum delay spread 
of each user k is much smaller than a symbol duration Tb. 

As mentioned above, during the detection phase 40, the 
number of degrees of freedom of the second signal y(t) is equal to the 
number K of symbols bk sent by the users k. The rate of innovation p of the 
15 second signal y(t) is then equal to its information rate K/Tb. The sampling 
frequency fs must then be 



which is still well below the chip rate IH'c of the received signal. 

In the case of noisy transmission, the propagation parameters ak^'^ 
20 and xk^'^ of the transmission channel c and consequently the transmitted 
symbols bk can not be estimated perfectly. In order to perform a better 
estimation of these channel parameters and therefore a less bit error rate, 
the sampling frequency fs should be higher than the rate of innovation p. 
The sampling frequency fs required for good estimation accuracy will then 
25 depend on signal-to-noise ratio (SNR), yet in most cases encountered in 
practice, that sampling frequency fs is still far below the chip rate 1/Tc. 

in a variant preferred embodiment of the method, the filtering 
signal (p(t) is a Gaussian signal. Since the Gaussian signal has an exponential 
decay, the problem of intersymbol interference, which can potentially arise 
30 with the use of a sine filtering signal, can be avoided. 
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In another variant preferred embodiment, the method according 
to the invention is applied to a CDMA system where the receiver r is 
connected to an array of antennas (figure 5). The array of antennas is 
composed of a regularly spaced aligned antennas i. After demodulation, 
5 the signal yj(t) received by each antenna i is a superposition of multiple 
copies of the signals transmitted by the K users k over the transmission 
channel c, each copy having a possibly different amplitude attenuation and 
a possibly different relative delays with respect to a reference at the 
receiver r. Each copy arrives at the antenna array at the incident angle 9k^'^ 
10 which is the same for all antennas i in the array. However, there will be a 
fixed phase difference of 

e ^ = e^*-^ 

between the copies received at each two consecutive antennas i, 
where D denotes the spacing between two adjacent antennas, ©c is the 
1 5 carrier angular frequency and c is the speed of light. This holds under the 
assumption that the carrier frequency is relatively high compared to the 
bandwidth of the transmitted signal, which is typically true in practical 
cases. 

The copies of the same signal component received by each two 
20 consecutive antennas i are thus shifted by a fixed phase delay that depends 
on the direction of arrival 9k^'^ of the specific signal component. Thus the 
same phase difference will exist between the Fourier series coefficients of 
the same signal component received by each two consecutive antennas i. 

During the training phase 30, each signal yj(t) received by each 
25 antenna i is, as in the previously described embodiments of the method, 
first filtered with a lowpass filter f and then sampled at a frequency fs 
higher than the innovation rate p of the signal yj(t), but lower than the 
frequency given by Shannon's sampling theorem. As previously described, 
sets of spectral values Yj[m] are then computed from each generated set of 
30 sampled values yilnTsL for computing a matrix Di with each said set of 

spectral values Yi[m], with the known values of the bits bk^^^ of the training 
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sequences bkt and with the known spectral values Sk[rr\] of the coding 
sequences Sk(t}. The elements dmn of the computed matrices D,- are 
therefore given by 



1=1 

where Ag denotes the antenna gain (assumed to be the same for 
all antennas). The unknown propagation parameters Tk^^9k^'^ and ak^'^ of the 
user k can then be estimated by considering only the k-th columns of the 
matrices Di. Thus, for each user k, a matrix Fk made up of the k-th column 
of each matrix Dj is defined by 



10 F,=(D,(:,k) D,(:,k) D3(:,k) ... D,(:,k)) = 



d\k ^Ik 



V^Mk ^Mk ^Mky 



From the matrix Fk it is then solved for the propagation 
parameters Tk^'^,Ok^'^ and ak^'^ of user k, for instance using 2-D subspace 
methods for harmonic retrieval. Therefore, by repeating this operation for 
all users k, the unknown propagation parameters tk^'^Gk^'^ and ak^'^ of all 
15 users k can be estimated. The problem of multiuser parameter estimation is 
in this case decomposed into a series of 2-D estimation problems. 

The one skilled in the art will recognize that it is also possible to 
estimate all the propagation parameters Xk^^Ok^'^ and ak^'^ of the user k by 
considering only the first two columns of the matrix Fk, that is, by 
20 considering the signals from only two antennas i, and using for instance 1- 
D subspace methods, which could further reduce the computational 
requirements. However, this last variant embodiment of the method 
typically leads to less reliable estimation of the propagation parameters 
than the previous one. 
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In the detection phase, beams of antenna arrays are oriented 
according to the previously estimated directions of arrival and detection is 
performed for instance with the known 2-D RAKE method but on a lowpass 
filtered version of the received signals yi(t) sampled at a sampling frequency 
5 fs higher than the innovation rate p of the received signals yi(t), but lower 
than the frequency given by Shannon's sampling theorem. 

In a further preferred embodiment, the inventive method is 
applied to a CDMA communication system using a specific set of coding 
sequences Sk(t) that are orthogonal to each other on their lowpass filtered 

10 version. In this preferred embodiment the training phase 30 is performed 
as described in the previous embodiments. The detection phase 40 is 
performed by first filtering the received signal y{t) with a filter f and then 
sampling it at a frequency fs lower than the frequency given by Shannon's 
sampling theorem, but higher than the rate of innovation p of the received 

1 5 signal y(t), for generating a set of sampled values y[nTs]. These sampled 
values y[nT5] are then processed by a bank of digital matched filters, each 
matched to the sampled lowpass filtered version of a user's coding 
sequence Sk(t). As the lowpass filtered versions of the coding sequence Sk(t) 
of the users k are orthogonal to each other, the output of each matched 

20 filter determines the value of the information sent by the corresponding 
user k. 

The method according to the invention was demonstrated above 
for the case of a CDMA system with K users k transmitting over a multipath 
fading transmission channel c. The one skilled in the art will however 
25 recognize that the inventive method can also be applied to communication 
systems with transmission channels having a single propagation path as 
well as for single-user channel estimation. 

In the preferred embodiments described above, the method 
according to the invention and the related receiver using said method are 
30 applied to CDMA communication systems. The one skilled in the art will 
however recognize that the method according to the invention can be 
applied to other bandwidth-expanding communication systems, such as 
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ultra wideband systems (UWB). More precisely, UWB system can be roughly 
considered as a special case of CDMA systems where the signature sequence 
assigned to a user fs a very short pulse, much shorter than the bit duration 
Tb. In the typical practical case, this pulse is a first derivative Gaussian 
5 function. 

The typical implementation of an UWB communication system is 
a pulse position modulation system (PPM system) using a wideband pulse 
shape. In such systems, in order to recover the information contained in the 
received signal, it is necessary to estimate the locations of the transmitted 
10 pulses. The one skilled in the art will recognize that this task can be viewed 
as a particular case of the time-delay estimation described above during the 
training phase in the case of a CDMA system. 

During the training phase, a user sends a stream of regularly 
spaced pulses in order to allow for the receiver to estimate, in a similar 

15 manner to what is described above, the propagation parameters, that is, 
the relative time delays induced by the channel t^'^ as well as the 
propagation coefficients a^^\ During the detection phase, the receiver has to 
solve the same time-delay estimation problem as in the preceding training 
phase. In order to recover the information sent, the previously estimated 

20 time delays induced by the channel must be subtracted form the time 
delays estimated in the detection phase. 

The inventive method, in its various preferred embodiments, is 
described above for performing both channel estimation during the 
training phase and signal detection during the detection phase. It is 

25 however possible to implement the inventive method in a bandwidth- 
expanding communication system for performing only channel estimation 
and to use a method known from the prior art for performing signal 
detection, or to implement the inventive method in a bandwidth- 
expanding communication system for performing only signal detection and 

30 to use a method known from the prior art for performing channel 
estimation. 



